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Gravitational lensing is a powerful tool for the study of the distribution of 
dark matter in the Universe. The cold-dark-matter model of the formation of 



1 



large-scale structures predicts^^^ the existence of quasars gravitationally lensed 
by concentrations of dark matter^ so massive that the quasar images would 
be split by over 7 arcsec. Numerous searches*^"^^ for large-separation lensed 
quasars have, however, been unsuccessful. All of the roughly 70 lensed quasars 
known^^, including the first lensed quasar discovered^'^, have smaller separations 
that can be explained in terms of galaxy-scale concentrations of baryonic mat- 
ter. Although gravitationally lensed galaxies^^ with large separations are known, 
quasars are more useful cosmological probes because of the simplicity of the 
resulting lens systems. Here we report the discovery of a lensed quasar, SDSS 
J1004+4112, which has a maximum separation between the components of 14.62 
arcsec. Such a large separation means that the lensing object must be dominated 
by dark matter. Our results are fully consistent with theoretical expectations^"^ 
based on the cold-dark-matter model. 

For bright quasars (i-band magnitude i < 19), the probability of gravitational lensing^^ 
is only about 0.1%; the majority of these lenses have small separations, due to a single massive 
galaxy. The fraction of large-separation lensed quasars is predicted to be 0.01% or less^~^; 
thus it is not surprising that none have been found to date^~^^. In order to find such objects, 
we need samples of tens of thousands of quasars, such as generated by the Sloan Digital Sky 
Survey^^'^^ (SDSS). The SDSS is conducting both a photometric survey^*^"^^ using five broad 
optical bands^^ {u, g, r, i and z) and a spectroscopic survey^^ of 10,000 square degrees of the 
sky centered approximately on the North Galactic Pole, using a dedicated wide-field 2.5-m 
telescope at the Apache Point Observatory. 

We searched for large-separation lensed quasars in a sample of ~29,500 spectroscopically- 
confirmed SDSS quasars^^ at redshifts z of 0.6 — 2.3, a sample roughly three times larger 
than those used in previous searches. Even with this large sample, the expected number of 
large-separation lensed quasars is of the order of unity. In the field around each quasar in 
the sample, we searched for stellar objects with colours differing by less than 0.1 from those 
of the quasar, with separations between 7'.'0 and 60'.'0 and with flux greater than one-tenth 
that of the quasar. SDSS J1004-I-4112 was identified as a 'quadruple' large-separation lensed 
quasar candidate using these criteria. Only one of the four components (component B, see 
below) has an SDSS spectrum (the SDSS hardware^^ does not allow pairs of objects sepa- 
rated by less than 55" to be observed on a single plate), and therefore, we obtained spectra of 
all four components using the Keck I telescope at the W. M. Keck Observatory. The results 
are shown in Fig. 1. All four components indeed show quasar-like features, with all emission 
lines giving a consistent redshift z = 1.734±0.002; the velocity differences of the quasar 
components are ~ 100 km s~^, comparable to the observational uncertainty. Although it 
is not obvious from Fig. 1, there are C IV absorption line systems at 2; = 1.732 in each 
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of the four quasar spectra; this is an absorption system associated with the quasar itself, 
further supporting the lensing hypothesis: the four quasar images are from the same physical 
source. The differences in their spectra may be explained by the modest time-variability of 
the source quasar over ~1 year, the expected gravitational lensing time-delay^^ among those 
different images. 

Additional strong support for the lensing hypothesis comes from the identification of 
the galaxy cluster responsible for the large-separation lensing. From the observed image 
separations (the maximum separation is 14'.'62, between images B and C), we infer that the 
lensing object should have a velocity dispersion in excess of 600 km s~^. Thus the lensing 
object cannot be a single galaxy, but must be rather a group or cluster of galaxies that 
has a sufficiently concentrated distribution of dark matter. To identify the lensing object, 
we obtained deep optical images of the system using the Subaru telescope of the National 
Astronomical Observatory, Japan. The result is shown in Fig. 2. A number of galaxies are 
clearly detected around component G, suggesting that it is the most luminous galaxy of the 
cluster. We obtained a spectrum of component G using the Keck I telescope. The spectrum 
shows a number of absorption features characteristic of a early-type galaxy at a redshift of 
z = 0.6799 ± 0.0001. We also obtained spectra of two faint galaxies immediately to the 
south-west of component G (Fig. 2) using the Faint Object Camera and Spectrograph^^ 
of the Subaru telescope. The redshifts of these two faint galaxies are z = 0.6751 ± 0.0001, 
strongly suggesting a cluster of galaxies at 2~0.68 centred on component G. Clusters are 
dominated by elliptical galaxies, which all have very similar spectral energy distributions. 

Many of the faint galaxies in Fig. 2 (~40 galaxies around component G) have colours 
that are similar to that of component G. The colours are consistent with the expected colours 
of elliptical galaxies at z~0.68 {g — r~1.8 and r — i~l.l). In addition, there is an X-ray 
source in this direction detected by the ROSAT All-Sky Survey^'' (0.236 counts per second 
in a 473-s exposure). The emission, however, comes most probably from the quasar, because 
the detected X-ray flux is too strong for typical clusters of galaxies at z = 0.68. Finally, we 
note two possible arclets (highly distorted images of background galaxies due to gravitational 
lensing) in Fig. 2 (marked as 'arc?') close to component D. If future observations confirm 
that the arclets are indeed lensed background galaxies, they will provide strong additional 
constraints on the total mass distribution of the lensing cluster. 

The lensing interpretation is further supported by a theoretical model of SDSS J1004+4112. 
We fitted the positions of the four quasar components with a singular isothermal ellipsoid 
(SIE) plus external shear model using lens modeling software^^. The best-fit model is il- 
lustrated in Fig. 3. The positions and relative brightnesses of all components agree well 
with the lens model predictions. The centre of the lensing mass is offset from the centre of 



3 



component G by about 10 kpc at the cluster redshift, but brightest galaxy of a cluster is not 
always found exactly at the centre of the potential well of that cluster. 

The identical redshifts {z = 1.734) and the spectral energy distributions of the four 
lensed components, the existence of a lensing cluster of galaxies {z = 0.68), and the presence 
of possible arclets confirm the hypothesis that the quasar is lensed by this cluster. Further- 
more, a theoretical lensing model involving the cluster and external shear simultaneously 
accounts for the observed geometry of the system and the relative brightness of the images. 
The present work represents the discovery of a long-predicted but previously undetected 
population of large-separation lensed quasars. 

1. Narayan, R. & White, S. D. M. Gravitational lensing in a eold dark matter universe. Mon. Not. R. 
Astron. Soc. 231, 97-103 (1988). 

2. Wambsganss, J., Cen, R., Ostriker, J. P. & Turner, E. L. Testing Cosmogonie Models with Gravita- 
tional Lensing. Science 268, 274-276 (1995). 

3. Keeton, C. R. & Madau, P. Lensing Constraints on the Cores of Massive Dark Matter Halos. Astro- 
phys. J. 549, L25-L28 (2001). 

4. Wyithe, J. S. B., Turner, E. L. & Spergel, D. N. Gravitational Lens Statistics for Generalized NEW 
Profiles: Parameter Degeneracy and Implications for Self-Interacting Cold Dark Matter. Astrophys. 
J. 555, 504-523 (2001). 

5. Takahashi, R. & Chiba, T. Gravitational Lens Statistics and the Density Profile of Dark Halos. 
Astrophys. J. 563, 489-496 (2001). 

6. Oguri, M. Constraints on the Baryonic Compression and Implications for the Fraction of Dark Halo 
Lenses. Astrophys. J. 580, 2-11 (2002). 

7. Navarro, J. F., Frenk, C. S. & White, S. D. M. A Universal Density Profile from Hierarchical Clus- 
tering. Astrophys. J. 490, 493-508 (1997). 

8. Maoz, D., Rix, H., Gal- Yam, A. & Gould, A. Survey for Large-Image Separation Lensed Quasars. 
Astrophys. J. 486, 75-84 (1997). 

9. Ofek, E. O., Maoz, D., Prada, F., Kolatt, T. & Rix, H. A survey for large-separation lensed FIRST 
quasars. Mon. Not. R. Astron. Soc. 324, 463-472 (2001). 

10. Phillips, P. M. et al. The JVAS/CLASS search for 6-arcsec to 15-aresec image separation lensing. 
Mon. Not. R. Astron. Soc. 328, 1001-1015 (2001). 

11. Zhdanov, V. I. & Surdej, J. Quasar pairs with arcminute angular separations. Astron. & Astrophys. 
372, 1-7 (2001). 

12. Kochanek, C. S. et al. CASTLES Survey, at <http://cfa-www.harvard.edu/castlcs/> (2003). 

13. Walsh, D., Carswell, R. F. & Weymann, R. J. 0957 + 561 A, B - Twin quasistellar objects or 
gravitational lens. Nature 279, 381-384 (1979). 



4 



14. CoUcy, W. N., Tyson, J. A. & Turner, E. L. Unlcnsing Multiple Arcs in 0024+1654: Reconstruction 
of the Source Image. Astrophys. J. 461, L83-L86 (1996). 

15. Turner, E. L., Ostriker, J. P. & Gott, J. R., Ill The statistics of gravitational lenses - The distributions 
of image angular separations and lens redshifts. Astrophys. J. 284, 1-22 (1984). 

16. York, D. G. et al. The Sloan Digital Sky Survey: Technical Summary Astion. J. 120, 1579-1587 
(2000). 

17. Stoughton, C. et al. Sloan Digital Sky Survey: Early Data Release. Astron. J. 123, 485-548 (2002). 

18. Gunn, J. E. et al. The Sloan Digital Sky Survey Photometric Camera. Astron. J. 116, 3040-3081 
(1998). 

19. Pier J. R. et al. Astrometric Calibration of the Sloan Digital Sky Survey. Astron. J. 125, 1559-1579 
(2003). 

20. Hogg, D. W., Finkbeiner, D. P., Schlegel, D. J. & Gunn, J. E. A Photometricity and Extinction 
Monitor at the Apache Point Observatory. Astron. J. 122, 2129-2138 (2001). 

21. Smith, J. A. et al. The u'g'r'i'z' Standard-Star System. Astron. J. 123, 2121-2144 (2002). 

22. Fukugita, M. et al. The Sloan Digital Sky Survey Photometric System. Astron. J. Ill, 1748-1756 
(1996). 

23. Blanton, M. R. et al. An Efficient Targeting Strategy for Multiobject Spectrograph Surveys: the 
Sloan Digital Sky Survey "TiHng" Algorithm. Astron. J. 125, 2276-2286 (2003). 

24. Richards, G. T. et al. Spectroscopic Target Selection in the Sloan Digital Sky Survey: The Quasar 
Sample. Astron. J. 123, 2945-2975 (2002). 

25. Oguri, M., Taruya, A., Suto, Y. & Turner, E. L. Strong Gravitational Lensing Time Delay Statistics 
and the Density Profile of Dark Halos. Astrophys. J. 568, 488-499 (2002). 

26. Kashikawa, N. et al. FOCAS: The Faint Object Camera and Spectrograph for the Subaru Telescope. 
Piibhcations of the Astronomical Society of Japan 54, 819-832 (2002). 

27. Cao, L., Wei, J.-Y. & Hu, J.-Y. High X-ray-to-optical flux ratio RASS-BSC sources. I. The optical 
identification. Astron. Sz Astrophys. Supplement 135, 243-253 (1999). 

28. Keeton, C. R. Computational Methods for Gravitational Lensing. Preprint astro-ph/0 102340 at 
<http://xxx.lanl.gov> (2001). 

29. Oke, J. B. et al. The Keck Low-Resolution Imaging Spectrometer. Piibhcations of the Astronomical 
Society of the Pacific 107, 375-385 (1995). 

30. Miyazaki, S. et al. Subaru Prime Focus Camera - Suprime-Cam. Publications of the Astronomical 
Society of Japan 54, 833-853 (2002). 



Acknowledgments Funding for the creation and distribution of the SDSS Archive has 
been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the Na- 
tional Aeronautics and Space Administration, the National Science Foundation, the U.S. 



5 



Department of Energy, the Japanese Monbukagakusho, and the Max Planck Society. The 
SDSS Web site is http:/ /www. sdss.org/. 

The SDSS is managed by the Astrophysical Research Consortium (ARC) for the Partic- 
ipating Institutions. The Participating Institutions are The University of Chicago, Fermilab, 
the Institute for Advanced Study, the Japan Participation Group, The Johns Hopkins Uni- 
versity, Los Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), 
the Max-Planck- Institute for Astrophysics (MPA), New Mexico State University, University 
of Pittsburgh, Princeton University, the United States Naval Observatory, and the University 
of Washington. 

This paper is based in part on data collected at the Subaru Telescope, which is oper- 
ated by the National Astronomical Observatory of Japan, W. M. Keck Observatory, which 
is operated as a scientific partnership among the California Institute of Technology, the 
University of California, and the National Aeronautics and Space Administration, and the 
Apache Point Observatory (APO) 3.5-meter telescope, which is owned and operated by the 
Astrophysical Research Consortium. 

Part of this work was performed under the auspices of the U.S. Department of Energy 
at the University of California Lawrence Livermore National Laboratory. 

Competing interests statement The authors declare that they have no competing finan- 
cial interests. 

Correspondence and requests for materials should be addressed to N.I. 
(e-mail: inada@utap.phys.s.u-tokyo. ac.jp) 



6 



Table 1. ASTROMETRY AND PHOTOMETRY FOR SDSS J1004+4112 



Object 


R.A.(J2000)^ 


Dec.(J2000)^ 












A 


10 04 34.794 


+41 12 39.29 


18.67+0.03 


18.70+0.02 


18.46+0.02 


18.43+0.05 


3'.'73 


B 


10 04 34.910 


+41 12 42.79 


19.05+0.06 


19.10+0.06 


18.86+0.06 


18.92+0.06 


O'.'OO 


C 


10 04 33.823 


+41 12 34.82 


19.71+0.03 


19.73+0.02 


19.36+0.03 


19.31+0.07 


14'.'62 


D 


10 04 34.056 


+41 12 48.95 


20.67+0.04 


20.51+0.04 


20.05+0.04 


20.00+0.13 


11'.'44 


G 


10 04 34.170 


+41 12 43.66 


22.11+0.40 


20.51+0.13 


19.54+0.09 


19.04+0.21 


8'.'44 



'^R.A., right ascension; Dec, declination. RA and Dec. are given in the J2000 system. These celestial 
coordinates were measured on the basis of the celestial coordinates of component B. The positional 
errors of components A, C, and D (not including the absolute positional errors of component B) are 
O'.'Ol and that of component G is 0'.'05 per coordinate. 

^g, r, i and z mean the magnitudes of each band. 

^Separation angles relative to component B 
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Figure 1 The Keck spectra of the four quasar components A-D and the brightest galaxy G in the 
lensing cluster. See Fig. 2 for these identifications (A-D, and G). The data were taken using the 
Low-Resolution Imaging Spectrometer^^ (LRIS) of the Keck I telescope. The exposure times were 
900 s for each component. The dispersion is l.OOA pixel"^. The data were reduced in a standard 
method using IRAF (IRAF is the image reduction and analysis facility, distributed by the National 
Optical Astronomy Observatories, which are operated by the Association of Universities for Re- 
search in Astronomy, Inc., under cooperative agreement with the National Science Foundation). 
The black solid line, the red solid line, the green solid line, and the blue solid line represent the 
spectra of components A, B, C, and D, respectively. The vertical gray dotted lines (3323.6 A, 3818.7 
A, 4235.1 A, 5218.5 A, and 7651. 8A) represent the positions of emission lines of the respective ions 
redshifted to z = 1.734 of Lya (1215.67 A), Si IV (1396.76 A), C IV (1549.06 A), C III] (1908.73 A), 
and Mg II (2798.75 A), respectively. All emission lines are clearly at the same redshift. The orange 
solid line represents the Keck spectrum of component G at the same dispersion. The exposure time 
was also 900 s for component G. The vertical thinner gray dotted lines (3933.7 A, 3968.5 A, 4304.4 
A, and 5175.3 A) represent the positions of absorption lines of the respective ions redshifted to 
z = 0.680 of Ca II H&K (3933.7 A and 3968.5 A), G-band (4304.4 A), and Mg I b-band (5175.3 
A), respectively. There are no data below ~ 4900A in the spectrum of component G. 
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Figure 2 The gri composite Subaru image of the field around SDSS J1004+4112. The data were 
taken using the Subaru Prime Focus Camera'^'' of the Subaru telescope. The magnitude limit is 
i ~ 26.0. The central 60" square is shown in an expanded view. The four quasar components are 
marked as A, B, C and D, and the bright galaxy located between the four quasar components is 
marked as G. The separation between components A and D is 12'.'77, and that between components 
B and C is 14'.'62. The positions (J2000) and the magnitudes of the components A-D and the 
brightest galaxy (component G) between the four quasar components are summarized in Table 1. 
Many faint galaxies can be seen — their positions and colours are consistent with being members 
of a cluster {z = 0.68) centred on component G. Two possible arclets (marked as 'arc?') can also 
be seen. The seeing had a full-width at half-maximum of O'.'O. 
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Figure 3 The best fit lens model prediction compared with the observation. We used a lensing 
model of a singular isothermal ellipsoid (SIE) with external shear. The best-fit model has an 
Einstein radius of the SIE model ctg = 6'.'906 (corresponding to a velocity dispersion ~ 700 km sec~^ 
at the cluster redshift), magnitude and position angle of the shear 7 = 0.250 and 9^ = —60.925° 
(measured East of North), and ellipticity and its position angle e = 0.498 and 9e = 21.434° 
(measured East of North), with a source quasar position (AR.A., ADec.)=(— 7'.'124, — 0'.'574) and 
a centre of lensing mass (AR.A., ADec.)=(— 7'.'387, — 0'.'004) relative to the centre of component A. 
The black filled circles represent the observed positions of components A, B, C, D and G, and the 
red filled circles represent the predicted positions of components A-D. The green solid line is the 
position of the caustic in the source plane, and the blue dashed line represent the critical curve in 
the image plane. The small brown filled cross is the predicted position of the source quasar, and the 
small orange triangle is the predicted position of the center of the lens mass. The differences between 
the observed and modelled image positions are much smaller than the observational uncertainties. 
The flux ratios predicted from the model, B/A, C/A and D/A are 0.78, 0.43 and 0.22, respectively. 
The total magnification of the quasar images which is predicted by the model is 56.48. The 
predicted flux ratios are close to the observational results; B/A= 0.69±0.04, C/A= 0.46±0.02, 
and D/A= 0.25±0.01 (measured from the i band image). Microlensing by substructures and/or 
reddening by the Mg 11 absorption line systems that are seen in each spectrum might be the cause 
of the differences between the predicted flux ratios and the observations. 
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